GW approximation with LSDA+U method and applications to NiO, MnO, and V2O3 
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A GW approximation (GWA) method named U+GWA is proposed, where we can start GWA 
with more localized wave functions obtained by the local spin-density approximation (LSDA)-I-U 
method. Then GWA and U+GWA are applied to MnO, NiO, and V2O3 in antiferromagnetic phase. 
The band gaps and energy spectra show an excellent agreement with the experimentally observed 
results and are discussed in details. The calculated width of d-bands of V2O3 is much narrower 
than that of the observed one which may be a mixture of t|g multiplet and single electron t2g level. 
GWA or U+GWA does not work also in the paramagnetic phase of V2O3 and the reason for thois 
is clarified. The method of the unique choice of on-site Coulomb interaction is discussed in details. 
The criterion for whether we should adopt GWA or U+GWA is discussed and is assessed with the 
help of the off-diagonal elements of the self-energy. 

PACS numbers: 71.10.-w, 71.15.-m, 71.20.-b 
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I. INTRODUCTION 

Electron correlation is important for electronic proper- 
ties in wide variety of materials, where it causes drastic 
change in physical properties in strongly correlated elec- 
tron systems with small change in electron/hole doping 
or an external field. Standard electronic structure cal- 
culations arc now based on the density-functional theory 
(DFT) and have progressed greatly the understanding 
of electronic structures in condensed matters. The po- 
larization function in DFT is derived by the linear re- 
sponse theory of homogeneous/inhomogeneous electron 
gas, where the electron-electron correlation is calculated 
with the random-phase approximation (RPA) and effects 
of static screening are included. However, the results 
could not be free from self-interaction, especially in the 
local (spin-)density approximation (L(S)DA)fi On the 
other hand, GW approximation (GWA) does not suffer 
from the problem of the self-interaction and can include 
effects of dynamical screening in the framework of many- 
body perturbation theory;^ Moreover, one can estimate 
the screened Coulomb interaction from the first princi- 
ples. With a recent progress in computational algorithm 
and computer facilities, GWA can be applied now to re- 
alistic materials. 

GWA is the first term approximation of the many-body 
perturbation series and the self-energy is replaced by the 
lowest order term of the perturbation expansion. GWA 
for realistic condensed matters is formulated usually with 
the LSDA Hamiltonian as an unperturbed one with- 
out self-consistent calculation of the one-body Green's 
function G and the LSDA exchange-correlation potential 
is subtracted later from the resultant GWA self-energy. 
The GWA self-energy is then given as 



where Go is the unperturbed Green's function. W is 
the dynamically screened Coulomb interaction with the 
random phase approximation (RPA), which is written as 



W 
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where v is the bare Coulomb interaction and is 
the lowest order irreducible polarization function x'' = 
— iGoGq. 

One possible way is that the quasi-particlc energy is 
determined self-consistently in the lowest GWA equation 
with which is called the "eigenvalue-only (e-only) self- 
consistency" ^ In fact, the e-only self-consistent approach 
may be good enough if LSDA would give reasonably good 
starting wave functions, though it may not be always the 
case. 

The self-consistent calculation of G in GWA may be 
another approach, but that without the vertex correc- 
tion leads to unphysical structure of spectra, e.g., too 
wide band width and disappearance of plasmon satel- 
lite, and violation of the /-sum rule, though it ensures 
conservation of particle number and energy and then an 
accurate total energy4i^ Other trials of partially self- 
consistent treatment of GWA, named "the quasiparticle 
self-consistent" GWA,— have been proposed to improve 
the quasiparticle band structure. The essence of these 
methodologies is in how to obtain localized wave func- 
tions in transition-metal oxides. 

Another possibility would be the establishment of a 
methodology to start with some unperturbed Hamil- 
tonian which gives localized wave functions or correct 
eigen-cnergies. In the present work, we propose a 
methodology of GWA, starting with wave functions of 
the LSDA+U mcthod^i^ in order to have localized wave 
functions and apply it to several transition-metal oxides, 
MnO, NiO and ¥203.^^ In Sec. HH the present method- 
ology, named U+GWA, is explained briefly. Section IIIII 



2 



is devoted to the application to antiferromagnetic insu- 
lators (AFI) MnO and NiO. Two systems are very typi- 
cal: The electronic structure of MnO can be obtained by 
GWA and, on the contrary, that of NiO by U+GWA. The 
physical reason for this choice is discussed. The method 
of unique choice of the on-site Coulomb interaction U is 
analyzed in details. We show that the starting wave func- 
tions by LSDA are satisfactory in MnO and are largely 
improved in NiO by LSDA-I-U with a proper value of U. 
The present U-hGWA is then appHed to AFI V2O3 in 
Sec. IIVI Summary and conclusion are given in Sec. |Vl 

II. U+GWA: GW APPROXIMATION 
STARTING FROM LSDA+U 



A. U+GWA 



where Ckn is the LSDA+U eigen-energy. The self-energy 
correction in Eq. ([5]) can be written as ASknl^-kn) = 
(i/'kn|S(.Ek„) - Vl'^sda - ^LSBA+v iV'kn)- We actually car- 
ried out energy-only self-consistent calculation to satisfy 
Eq. ®. 

Usually GWA causes an appreciable mixing between 
the LSDA valence and conduction bands, if we include 
the off-diagonal elements of the self-cnergyiii In the 
present work, the Coulomb U in the LSDA+U calculation 
is introduced in order to obtain localized wave functions. 
The choice of U will be discussed further in Sec. IIII Dl 
Our criterion for the value U in U+GWA is that the 
off-diagonal elements of the self-energy become negligi- 
bly small, which will be discussed in Sec. IIII El The 
exchange interaction J is set to be the value of the con- 
strained LSDA calculationii^ If the calculation of Green's 
function were done self-consistently, the results would not 
depend on the choice of values of U and ,/. 



We will present here a theoretical method of GWA 
with wave functions obtained by the LSDA+U method^ 
and we apply this method to NiO, MnO, and V2O3 in 
the proceeding sections. Present calculation is based 
on the linear muffin-tin orbital method with the atomic 
sphere approximation (LMTO-ASA) (Ref. [13) and also 
the LSDA+U method with Coulomb interactions of ro- 
tational invarianccji The LSDA+U method has shown 
a reasonable success for many d- and f-electron systems 
in the broad feature of electronic structure. Since the 
LSDA+U method is a kind of the static limit of GWA, it 
could be a good starting approximation of GWA4 The 
energy gap becomes larger in LSDA+U and the polariza- 
tion becomes smaller because of the energy denominator 
in the RPA polarization function. Moreover, wave func- 
tions of the LSDA+U method are more localized than 
those by LSDA because of strong on-site Coulomb inter- 
action and, therefore, can be more preferable for materi- 
als with strong electron-electron correlation. 

Once we start from the LSDA+U Hamiltonian, the 
exchange-correlation potential and the Hubbard terms 
of Coulomb interaction should be subtracted from the 
GWA self-energy as 

= ^ ^ ^LSDA ^ ^LSDA+U' (3) 

where Uj^gDA ^^1*^ exchange-correlation potential and 
Ulsda-i-u is the potential correction derived from the 
Hubbard term in the LSDA+U methodii This we call 
LSDA+U+GWA or U+GWA. In the present formula- 
tion, once we set U = and J ~ 0, the last term in 
Eq. ^ vanishes and U+GWA is reduced to GWA. 
Green's function G is defined as 

G{E) = iE-Ho- AI](^))-i (4) 

and the quasiparticle energy i5kn should be calculated by 
the self-consistent equation (the e-only self-consistency) 

Ewn — Ckri + R.cAI]kn(i?kn), (5) 



B. Coulomb interactions 

The effective Coulomb and exchange integrals between 
orbitals (pa and U{a,b) and J[a,b), are represented 
by the Racah parameters A, B, and C or by the Slater 
integrals Fq, F2, and F4. The Coulomb and exchange 
interactions within and between tig and Cg orbitals are 
as follows in the cubic symmetry: 

«t., = c/(e,e) = t/(rM?)-c/(c,c) 

= ^ U{u,u) = U{v,v) 

= A + 45 + 3C Fo + 4F2 + 36F4, (6) 

= A-2B + C = Fq-2F2'~4F4, (7) 
<^ = U{u, v)^A-4B + C = Fo~iF2+ 6F4, (8) 

u" = - V y U(a, b)^A + C^Fo- I4F4, (9) 
6 ^ — ' ^ — ' 

Jt2, = JiC,v)^J{vX)^JiC,0 

= 3B + C = 3^2 + 2OF4, (10) 
jcg = J{u,v) = 4B + C ^4F2 + 15F4, (11) 

= 2B + C = 2F2 + 25^4, (12) 

We have an experimentally observed relationship^ 
Fi/F'i — 0.63/9 = 0.07 and the above expressions are 
rewritten in terms of the Coulomb and exchange param- 
eters U and J with help of the relationships U = Fq and 
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TABLE I: The input Coulomb and exchange interaction parameters, U (eV) and J (eV), static limit of screened Coulomb 
interaction W{0) (eV), direct band gap i?G:d (eV), indirect band gap iJciid (eV), and the spin magnetic moment M{fiB) for 
MnO and NiO. The calculated direct and indirect band gaps are estimated from the calculated quasiparticle energy bands. 
The spherical average values of the bare Coulomb interaction (v) are 23.6 eV in MnO and 27.9 eV in NiO. We have found two 
U+GWA solutions in NiO with U — 1 eV and 2.0 eV. We believe that the solutions connecting to the ones of [/ > 2.5 eV 
continuously are correct. See the text in Sub. IIII Dl 
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The value of VF(0) depends on the orbital components and its averaged one is shown here. 
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J= 1(F2 + 9F4) as 



"t2g 


= Ue, =U+1.UJ, 


(13) 




= U-OAOJ, 


(14) 


1 


^ J7-0.63J, 


(15) 


3t2g 


= 0.77J, 


(16) 


JCg 


= 0.89J, 


(17) 


u" 


= C/-0.17J, 


(18) 


.]" 


= 0.66J. 


(19) 



These U and J may be the input parameters of LSDA+U 
method or the output Coulomb and exchange interac- 
tions in U+GWA method. 

The energy difference in one-electron spectra between 
an electron affinity level Ta and an electron ionization 
level Ti is evaluated as 

{ETiTa) - Et{G)} - {Et{G) - Erim 
^ Eg + {EciTa) - EciG)} - {EciG) - EdT,)}, 

(20) 

where the total and Coulomb energies of a multiplet T 
are written as Et{T) and Ec(r) respectively, Et{G) is 
the total energy of the ground state and Eq is the one- 
electron energy gap without the Coulomb interactions, 
e.g. a hybridization gap etc. In order to evaluate the 
Coulomb and exchange interactions, we will assume, in 



later sections, the form 

i [U{m,m')n''„^n;^ 

+{(7(m,m') - J(m,m')}«„,], (21) 

where nj^^ is the electron occupation. 

We will analyze the calculated energy spectra by using 
these formulas and show that the Eq vanishes in all these 
materials MnO, NiO and also V2O3. 



III. ELECTRONIC STRUCTURE OF NIO AND 
MNO BY GWA AND U+GWA 



A. Details of calculations 

We use 4x4x2 k-point mesh in the Brillouin zone of 
MnO and NiO. The set of the maximum orbital angular 
momentum of the LMTO basis in Mn, Ni, O and empty 
spheres are chosen to be (ffdp) . In the calculation of the 
self-energy, the product basis scheme is used^^ii^ and the 
maximum total angular momentum of the product bases 
is set to be (ffdp) for the calculation of the correlation 
part S'^ of the self-energy. This choice reduces the num- 
ber of the product bases used in S"^ from 2,376 to 360. 
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FIG. 1: (Color online) (a) Quasiparticle energy bands along 
symmetry lines (F = (0,0,0), X=^(1,0,0), K=2i(|,|,0) 
and point 22L(i^ q) equivalent to X-point) and (b) quasipar- 
ticle density of states of antiferromagnetic MnO. Solid and 
broken lines refer to those by GWA and those by LSDA. The 
energy zero-th is fixed at the top of the valence band. 



The e-only self-consistency is achieved with the iterative 
procedure of five or ten times. 



B. Electronic structure of antiferromagnetic MnO 
by GWA 

The LSDA band gap of antiferromagnetic MnO is 
0.80 eV in contrast with the observed one of 3.6~3.8 eV. 
The observed energy interval between the principal 
peaks of x-ray photoemission spectroscopy (XPS) and 
bremsstrahlung isochromat spectroscopy (BIS) is approx- 
imately 9.5 cV and corresponds to 



i?T((4j^(eJ)2)} 



{i?T((4,)^(eJ)^t^^))-^xu^2,y 
= ^G:MnO + U + 2.57 J, (22) 




Energy [eV] 



FIG. 2: (Color online) Imaginary part of the Green's functions 
(l/7r)|ImG(i;)| of antiferromagnetic MnO by GWA. (a) The 
total Green's function (solid curve) and experimental XPS 
and BIS spectra (dotted line)^^. (b) The partial Green's func- 
tion per atom. The energy zero-th is fixed at the top of the 
valence band. 



where i?G:MnO is the hybridization gap in MnO. 

The spin polarization in d-orbital is almost complete 
and the calculated numbers of d-electrons of majority and 
minority spins are 4.73 and 0.42, respectively. Since the 
d-d transition for the RPA polarization should be small 
because of almost filled d-band of the majority spin and 
empty d-band of the minority spin. Once we adopt GWA 
for MnO, the screening effect on the Coulomb interaction 
is small and the Coulomb interaction cannot be largely 
screened. Therefore, it is expected that the energy gap 
should be open appreciably in MnO with GWA. 

In Fig. [2 the quasiparticle band structure by GWA 
with the e-only self-consistency is depicted along the [110] 
direction from F-point (E- and S-lines) and [100] direc- 
tion (A-line). It must be noticed that the point ^(1, 1, 0) 
along the extended direction of S-line (1,1,0) is equiva- 
lent to the X-point ^(1,0,0). The corresponding quasi- 
particle density of states (DOS) is shown in Fig. [1] (b), 
compared with DOS by LSDA. 

In fact, GWA starting from LSDA gives rise to the 
band gaps to be 3.05 eV (indirect gap) and 3.51 eV (di- 
rect gap), in good agreement with the experimentally ob- 
served one Eq ~ 3.6 ^ 3.8 eV. The magnetic moment is 
evaluated to be M = 4.33 /xb comparable to the observed 
one M = 4.6 ~ 4.9 fiB- Therefore, we believe that our 
expected mechanism for less screened Coulomb interac- 
tion is the case in MnO. The U-I-GWA is also applied to 
AFI MnO with the c-only self-consistency as summarized 
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in Tabic U in the range of i7 = ~ 2 eV and J = 0.86 eV 
and the spectra and the quasiparticle band structure are 
much the same as those of GWA. We will discuss, in 
Sees. IIII Dl and UlI El the choice of U value and the start- 
ing wave functions in GWA or U+GWA of MnO. 

Figure [5] depicts the imaginary part of the Green's 
function (l/7r)|ImG'(£')| by GWA. The highest occupied 
band is the one with strongly hybridization between O p 
and Mn d (eg) states of majority spin and the lowest un- 
occupied one is mainly Mn d (t2g) states of minority spin. 
This fact tells that little d-d transition does contribute 
to the screening phenomena and also that MnO is just 
in the midway between of the charge-transfer type and 
Mott-Hubbard type insulators. The static limit of the 
screened Coulomb interaction W{0) is calculated to be 
7.07 eV, which means that the screened effect is small. 
We can conclude that GWA can reproduce the energy 
spectrum of AFI MnO satisfactorily, since little d-d tran- 
sition of the polarization can participate in the screening 
phenomena. 



C. Electronic structure of antiferromagnetic NiO 
by U+GWA 

NiO is one of materials where the LSDA wave func- 
tions may be extended in spaceJ^ The LSDA spectrum 
shows that the majority spin d-bands are almost full and 
the minority ones are half-full and that the band gap 
Eq is 0.11 eV (indirect gap) and 0.49 cV (direct gap) in 
contrast with the experimentally observed one of about 
4.0 ~ 4.3 eV. 

The observed energy interval between the principal 
peaks of XPS and BIS is approximately 6 ~ 6.2 eV and 
corresponds to 



{i?T((tig)^(t^g)^(cl )^e^)) - i?T((t^g)^(t^g)^(el )^)} 



G:NiO 



G:NiO 



U + 0.26 J, 



(23) 



where Eg-.NiO is the hybridization gap in NiO. 

GWA starting from LSDA gives rise to the band gap 
of Eg — 0.21 eV and the screened Coulomb interaction is 
estimated as W{0) = 1.51 eV, which does not agree with 
the experimental results. This is because the polarization 
function is too large due to small energy denominator 
of the RPA polarization function. We would expect a 
substantial change in wave functions and the polarization 
function led by an increase in the band gap, once we 
use GWA procedure under the condition of opening the 
gap.— The LSDA-I-U method gives better spectrum in 
NiO^^ 

U-I-GWA is applied to AFI NiO in the present work 
with the e-only self-consistency and a part of the re- 
sults is shown in Table HI compared with those by LSDA, 
LSDA+U, GWA and experiments. In fact, we have used 




Q -10 
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FIG. 3: (Color online) (a) Quasiparticle energy bands and (b) 
quasiparticle density of states of antiferromagnetic NiO. Sym- 
metry points and lines are the same as those in Fig. [T] Solid 
and broken lines refer to those by U-I-GWA and LSDA-I-U 
with U — 2.5 eV and J = 0.95 eV, respectively. In (a), the 
cross marks refer to the result by the angle resolved photoe- 



mission spectra.— 
the valence band. 



The energy zero-th is fixed at the top of 



two initial guesses of eigen-energics for the e-only self- 
consistent calculation for U ~ 1.0 eV and U = 2.0 cV; one 
starts from the LSDA-I-U eigen-energies and the other is 
from the LSDA+U eigen-energies plus some additional 
shift. Then we found different converged solutions for re- 
spective initial guesses. The cases U — and U > 2.5 cV 
give unique solutions in two starting eigen-values. The 
calculated values of these solutions arc summarized in 
Table [D We believe that the results connecting contin- 
uously to solutions oi U > 2.5 eV should be the correct 
ones and the other ones must be artifact. 

The calculation has been carried out in the range of 
values of C/ = - 8.0 cV and J ^ 0.95 eV (the value 
by the constrained LDA). The calculated peak position 
(at about 5 eV) of unoccupied bands with U = 2.5 eV is 
in an excellent agreement with experiments. The off- 
diagonal elements of the self-energy are extremely re- 
duced a.t U ~ 2.5 eV. Then, we adopt the case of 
U ~ 2.5 eV as appropriate choice. This choice will be 
discussed in more details in Sec. IIII Dl The band gaps 



6 



[l/eV] 
10 





Ni-1 



2g 



Ni-e, 





is.'! 



Energy [eV] 



10 



FIG. 4: (Color online) Imaginary part of the Green's functions 
(l/7r)|ImG(-B)| of antiferromagnetic NiO by U+GWA with 
U = 2.5 eV and J = 0.95 eV. (a) The total Green's function 
(solid line) and experimental XPS and BIS spectra (dotted 
line).— (b) The partial Green's function per atom. The energy 
zero-th is fixed at the top of the valence band. 



are 3.97 eV (direct) and 3.46 eV (indirect) and the mag- 
netic moment M = 1.46 /xb are in a good agreement with 
experimental ones as shown in Table [J The static limit 
of the screened Coulomb interaction W{0) is 6.03 eV. 

In Fig. [21 the quasiparticle band structure by U+GWA 
with U = 2.5 eV is depicted along the [110] and [100] 
symmetry lines, compared with the angle resolved pho- 
toemission spectroscopy data^^ and this shows an ex- 
cellent agreement between them. The d-band width 
in Fig. [3]^ a) (the energy width of upper occupied eight 
bands) becomes wider than that of LSDA (but does not 
change much from that of LSDA-f U), since the occupied 
O p-bands come upward and the hybridization mixing be- 
tween Ni d and O p becomes larger. It should be noted 
that the Ni 4s-band comes down slightly below the flat 
unoccupied Ni 3d-band. The corresponding quasiparti- 
cle DOS is shown in Fig. ^h), compared with DOS by 
LSDA+U. 

Figure [5] shows the imaginary part of Green's func- 
tion (l/7r)|lmG(£:)| by U4-GWA with U = 2.5 eV, with 
the photoemission and inverse photocmission spectral 
The most striking feature of i|ImG| is the band width 
and the intensity distribution in the spectra, e.g. the 
weight of the spectrum is shifted to the upper region of 
the occupied bands. The weight of O p orbitals in the 
highest occupied bands increases due to the stronger hy- 
bridization between O p and Ni d states and the lowest 
unoccupied ones arc mainly Ni d states. This fact is in 
a good agreement with observed one and NiO is of the 
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FIG. 5: (Color online) U dependence of the static screened 
Coulomb interaction Wu- x: LSDA-I-U, open square □: 
U-fGWA, broken line: Eq. (I25j and closed square: Eq. ^7} 
for MnO. +: LSDA-^U, open circle Q: V+GWA, solid line: 
Eq. (ESJ, and closed circle: Eq. ^ for NiO. See Subsec HITCl 
for the explanation on the two solutions of U+GWA (Q) at 
(7 = 1 and (7 = 2 in NiO. 



charge-transfer typci^ 

The screened Coulomb interaction can be calculated in 
GWA and U+GWA as shown in Table HI The screening 
effects are strong in the energy region < < 20 eV in 
GWA, since all d-electrons participate in the screening 
phenomena. The U+GWA can give much larger 1^(0) 
and reproduce the energy spectrum of antiferromagnetic 
NiO very satisfactorily. 



D. Unique choice of U in U+GWA in MnO and 
NiO 

We discuss, in this subsection, the U dependence of the 
static screened Coulomb interaction W{0) and explain a 
unique choice of U values in MnO and NiO. 

The irreducible polarization function in LSDA+U 
method may be expressed approximately in the static 
limit as 



, 0:LSDA-|-U 
A(7 



(24) 



where {E}P-^) is the effective band gap of the transition- 
metal d-bands of the minority spin by LSDA, and 6^ is a 
numerical factor of order 1 . The static screened Coulomb 
interaction Wu is then rewritten as 



W, 



LSDA+U 
U 



TX.LSDA. 



1 



1 



U 



(25) 



with the spherical average of matrix elements of the 
bare Coulomb interaction {v) and W^f^^ = {v)/{l — 
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XutT^^^i^))- Wc define the effective band gaps (E^P^) 
to be 3.7 eV for MnO and 1.4 eV for NiO from the partial 
density of states of the quasiparticle bands. We also es- 
timate the 6^ values as 0.665 for MnO and 0.764 for NiO 
so that W^^^ equals to the values of static screened 
Coulomb interaction by LSDA. 

The Wu values by LSDA+U and U+GWA methods 
are shown in Fig. [S] Equation (j25[) can represent very 
nicely the values of Wjj by LSDA+U in both MnO and 
NiO. In the MnO case, the values of LSDA+U method 
are shifted in the whole U range of U+GWA and con- 
tinuously vary. In NiO, on the contrary, there are two 
solutions at U = 1 eV and 2 eV, and then they vary 
smoothly. 

To understand the static screening Coulomb interac- 
tion of the U+GWA method we can establish 
a similar model of the irreducible polarization of Eq. 
as 



^0:U+GWA _ 

Xu — 



{Ec 



p^U+GWA • 



(26) 



where (Eq) is the hybridization gap appearing in 
Eqs. jSni), ^ and ^ and 6^ is the coefficient al- 
ready determined in Eq. (P^ . Then the static screened 
Coulomb interaction may be expressed as 



W, 



U+GWA 
U 



^0:U+GWA^^^ 



(27) 



The right-hand side of Eq. (gT]) is plotted in Fig. [5] by 
substituting the calculated static screened constant W{0) 
by U+GWA in the expression ;^^u+GWA ^-^^^ ^^^^ ^ 
and shows an excellent agreement between W{0) by 
U+GWA and Of course, though we cannot 

use the Eq. (|27|) to determine W{0) or input U value, 
this agreement implies that our treatment is consistent 
and, furthermore, the hybridization gap should be zero 
in both MnO and NiO. 

Though we cannot determine the value of U only 
in the framework of U+GWA, we can do with help 
of experimental spectra. We replace U in Eqs. ([H]) 
and ((23l) with ?7u+GWA = W{0) and we know now 
the values of the hybridization gap Eq equal to zero. 
With the help of the J value of the constrained LSDA 
calculation Ju-i-gwa, we have W(0) + 2.57Ju+gwa = 
W{0) + 2.21 ~ 9.5 i.e., W{0) = 7.29 eV for MnO and 
M^(0) + 0.26 Ju+GWA = W^(0) + 0.25 ~ 6.2 i.e., W{0) ^ 
5.95 cV for NiO. This is perfectly consistent with the re- 
sults obtained by the choices oi U = eV for MnO and 
U = 2.5 eV for NiO in U+GWA calculation. 



E. Legitimacy of U+GWA 

We have shown the calculated results of AFI MnO and 
NiO based on GWA and U+GWA, respectively. These 
two systems correspond to two different situations; 
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FIG. 6: (Color online) The value {V^:^ + 1 — a;)^ is shown, 
whose definition should be referred in the text and Eq. (|28|l . 
The vertical and horizontal axes are the LSDA eigen-energies 
for GWA or LSDA+U eigen-energies (U+GWA) ekn and ek„/ . 
The energy zeroth is set at the Fermi energy. The area having 
larger values in GWA disappears in U+GWA of NiO and no 
area having larger values exists in GWA of MnO. 



(1) The case (MnO) where the d-d transition contributing 
to the RPA polarization is small. The screening effect in 
GWA of the Coulomb interaction is relatively small and 
the screened Coulomb interaction is not reduced. In this 
case, large improvement can be expected in GWA calcu- 
lation. 

(2) The case (NiO) where the d-d transition contribut- 
ing to the RPA polarization is large. The screening ef- 
fect of the Coulomb interaction is large and the screened 
Coulomb interaction is reduced largely in GWA. In this 
case, large improvement could not be expected in GWA 
calculation and we should use U+GWA in order to make 
the energy gap open. 

Here we show how LSDA+U method improves the 
wave functions and polarization in NiO and how LSDA 
method provides reasonable results in MnO. We evaluate 
the mixing amplitude of the adopted wave functions by 
the off-diagonal element of the self-energy. With the off- 
diagonal matrix element of the self-energy ASknn', the 
two components of bands n and n' at k (the LSDA or 
LSDA+U eigen energies ekn and ekn') arc mixed with a 
ratio of (\A? + 1 — x)^ : 1,— where 



ekn + AI]kn(ekri) — Ekji' — AEkn^^kn') 



'({fkn + ekn'}/2) 



(28) 



When this ratio is small (cc is much larger than the unity), 
the resultant mixing is negligible. 

The values of this ratio for all pairs of bands ({kn} 
and {kn'}) are shown in Fig. [S) One can see, in NiO, 



that the area with larger ratio in GWA disappear in 
U+GWA of t/ = 2.5 eV. On the other hand, in MnO, 
no area with larger ratio exists in GWA. Therefore, in 
MnO, GWA does not promote a large mixing between or- 
bitals contributing to valence and conduction bands. On 
the contrary, in NiO, the U+GWA improves the start- 
ing wave functions and large off-diagonal elements of the 
sclf-encrgy disappear between orbitals contributing to va- 
lence and conduction bands. 



IV. ELECTRONIC STRUCTURE OF V2O3 BY 
U+GWA 

The antiferromagnetic V2O3 is the other example 
where a strong on-site Coulomb interaction plays a cru- 
cial role in the electronic structure. The crystal struc- 
ture in the low-temperature AFI phase is monoclinic and 
that in paramagnetic metallic (PM) phase above 150 K 
is corundum structure. The spin structure in AFI is de- 
picted in Fig. Efa), where each V'^+ ion has one spin- 
parallel (7 in Fig. [Tlja)) and two spin-antiparallel neigh- 
bors (/9) in the same layer perpendicular to the c-axis and 
one spin-parallel neighbor (a) in a different layer. We can 
have another view that V'^+ ions on one layer parallel to 
the c-axis are all ferromagnetically aligned and the inter- 
layer coupling is antiferromagnetic. The spin and orbital 
magnetic moments are observed to be (25") = 1.7 /is 
and {L) = —0.5 /is and the total spin is supposed to be 
S = 1^ The observed band gap is Eq = 0.35 ~ 0.66 eV. 
Under the trigonal symmetry around V ions, V dta level 
is split into non-degenerate &ig and doubly degenerate e^ 
levels. 

Mott pointed out that the metal-to-insulator (MI) 
transition in V2O3 at about 150 K is due to the electron- 
electron correlations 1^ The LSDA calculation shows that 
the aig states split into bonding and anti-bonding and 
that the lowest V d-state is aig bonding statci^ Then 
it was proposed that two aig electrons on a pair of V'^+ 
ions formed a spin-singlet and the remaining e^ electrons 
results a Mott-Hubbard model of an S" = i spin on each 
V'^+ ioui^ Ezhov et. al. studied its electronic structure 
by using the LSDA+U method with U = 2.8 eV and 
J = 0.93 eV, where the lowest V d-state is e^ and the 
aig is the next lowest. They showed that the orbital oc- 
cupation is predominantly e^Cg (a spin S = 1) with a 
small fraction of aigC^, and the band gap Eq — 0.6 eV 
in antiferromagnetic phase.— The estimated values for 
the on-site Coulomb and exchange interactions, U and 
J, are 2.8 and 0.93 eV, respectively. In the paramagnetic 
case, the electronic structures by LDA+U calculation is 
quite different from that in the antiferromagnetic case, 
e.g., metallic and large bonding-antibonding splitting of 
aig orbital which causes the majority configuration to be 
aigCg 1— 

LDA+DMFT (dynamical mean field theory) calcula- 
tion in paramagnetic phase using the Wannier-type or- 
bitals based on NMTO (N-th order muffin-tin orbital 



8 




FIG. 7: (Color online) (a) Spin structure of low temperature 
antiferromagnetic insulator phase of V2O3. Arrows stands for 
the spin-directions of V ions. The Greek letters a, /3, 7 and 
S denote the first, second, third and fourth nearest neighbor 
V-V pairs, respectively. The vectors a, b and their perpen- 
dicular vector c are the primitive vectors of the hexagonal 
unit cell and, the vectors a™, bm and Cm are those of the 
monoclinic cell of antiferromagnetic phase, (b) The Brillouin 
zone of the paramagnetic (solid lines) and antiferromagnetic 
(broken lines) phases of the corundum structure. The band 
structure in Fig. |S] is shown along the symmetry lines indi- 
cated by dot-dashed lines. The k-points L, Z and F are on 
the Brillouin zone of paramagnetic phase and F/2 on that of 
antiferromagnetic phase. 



method) formalism^^ gives very comprehensive physical 
picture of V2O3 with U = 4.2 eV, J = 0.7 eV and the 
U' ^ U — 2J ^ 2.8 eV, which presents a narrow reso- 
nance peak of a width of 0.5 eV and a broad lower Hub- 
bard bandit A crystal field splitting between aig and e^ 
is enhanced by the Coulomb interaction, then aig bands 
locates in the large Hubbard gap of e^ orbitals and the 
insulating gap is between the lower e^ Hubbard band and 
aig band. Therefore, the insulating gap in the paramag- 
netic phase is the crystal field gap but presumably not 
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TABLE 11: The Coulomb and exchange interactions, U (eV) 
and J (eV), static limit of screened Coulomb interaction M^(0) 
(eV), direct band gap -EG:d (eV), indirect band gap -Eciid 
(eV), and the spin magnetic moment M{fj,B) for V2O3. The 
calculated direct and indirect band gaps are estimated from 
the calculated quasiparticle energy. 





U 


J 


V2O3 

W^(0) EG:i 


Eg Ad 


M 


LSDA 








(metal) 


(metal) 


1.40 


LSDA+U 


2.6 


0.9 


2.62 


0.279 


0.109 


1.58 


GWA 








(metal) 


(metal) 


1.33 


U+GWA 


2.6 


0.9 


3.21 


0.943 


0.835 


1.93 


Constrained 


2.8 


0.9 










LSDA" 














exp. 








0.35, 


0.66' 


1.7" 



The value of VF(0) depends on the orbital components and its 
averaged one is shown here. 

"I.Solovyev et al. Phys.Rcv. B53, 7158 (1996). 

'G.A.Sawatzky and D.Post, Phys. Rev. B20, 1546 (1979); 
G.A.Thomas et al. Phys. Rev. Lett. 73, 1529 (1994). 

'^Rcfercnca^L 



the Hubbard gap. 



A. Results for antiferromagnetic V2O3 by U+GWA 

We use 4x4x4 k-point mcsli in the Brillouin zone of 
V2O3. The set of maximum angular momentum of the 
LMTO basis in V, O, and two different kinds of empty 
spheres are chosen to be (fdds). Those of product basis 
are (fdds), too, for the calculation of T,'^. The number 
of product bases used in E'^ is reduced from 10,728 to 
1,584. The two crystal structures, monoclinic and corun- 
dum structures, do not cause any significant difference 
in the electronic structures in LSDA and LSDA+U cal- 
culations for both paramagnetic and antiferromagnetic 
phases. Then, we use the corundum structure through- 
out the present work. 

The GWA results of V2O3, starting from LSDA, give 
rise to a metallic antiferromagnetic ground state, com- 
pletely different from the experimental situation. Both 
in LSDA and in LSDA+U, the d-bands of majority spin 
are partially filled and those of minority spin are empty. 
Though the screening effects in GWA must be large, the 
polarization function is too large due to small energy de- 
nominator, similarly to NiO. Therefore, we start calcu- 
lations with LSDA+U to prepare improved wave func- 
tions and polarization function and adopt U+GWA with 
the e-only self-consistency. The parameters U and J 
are set as 2.6 eV and 0.9 eV, respectively. The value 
of [/ = 2.55 eV may be the lower limit for opening the 
band gap in LSDA+U and one would get the metallic 
antiferromagnetic ground state for U < 2.55 eV. The 
parameters and calculated results are tabulated in Ta- 
ble [TTl The Brillouin zones of the paramagnetic and an- 
tiferromagnetic phases of corundum structure are shown 
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FIG. 8: (Color online) (a) Quasiparticle energy bands and 
(b) the density of states of antiferromagnetic V2O3. Solid 
and broken lines refer to those by U+GWA and LSDA+U 
with U = 2.6 eV and J = 0.9 eV. The energy zero-th is fixed 
at the top of the valence band. 



in Fig.EJb). 

Figure [SJa) depicts the quasiparticle band structure 
by U+GWA with U = 2.6 cV and J = 0.9 eV. The 
k-point F locates on the surface of the Brillouin zone of 
paramagnetic phase and the mid-point between F and F 
is denoted as F/2 locating on that of antiferromagnetic 
phase. The oxygen p-bands with a broad width of 3.7 eV 
locate at 6 eV below the Fermi energy which are not 
shown here. The d-band width becomes narrower due to 
the formation of a large energy gap. The corresponding 
quasiparticle DOS is shown in Fig. [5] (b), compared with 
DOS by LSDA+U. The e^ bands have a large exchange 
splitting, that of majority spin is almost fully occupied 
and that of minority spin is empty. Furthermore, the 
aig band of majority spin locates in the gap between e^ 
majority- and e^ minority-spin bands. 

The total Green's function i|ImG| is compared with 
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FIG. 9: (Color online) Imaginary part of the partial Green's 
functions of antiferromagnetic V2O3 by U+GWA with U — 
2.6 eV and J = 0.9 eV. (a) The total Green's function (solid 
line) and experimental XPS spectrum (dotted line).— (b) The 
partial Green's function per atom. The energy zero-th is fixed 
at the top of the valence band, aig and e^ are originated from 
dtig orbitals in the cubic symmetry and Cg is originated from 
dcg orbitals. 



the observed XPS spectra in Fig. [Hfa). The position and 
width of O p-bands are well reproduced in the calcu- 
lation but the width of the occupied V d-bands is too 
narrow. Figure [UJb) shows the imaginary part of the 
partial Green's function. The valence band just below 
the Fermi energy is mostly of c^ with small mixing of aig 
orbitals and the lowest conduction band is of ai^ char- 
acter. The ground-state configuration is, then, c^'^e^'^ 
of majority spins a and the e^"*^ band of minority spin 
— (T locates at higher energies by about 3.5 eV. The un- 
occupied aig band of majority spin is lifted up at the 
midway between majority- and minority-spin e^ bands. 
With an assumption of vanishing hybridization gap, the 
energy difference between main peaks of occupied e^ and 
unoccupied aig of parallel spins is estimated, replacing U 
with C/u+GWA = W^(0) and J by Ju+gwa, as 

-{i5c(efe^T)_i5c(e-T)} 

= - JtJ - - JtJ = W{0) - 1.17JU+GWA 

~ 2.16 eV, (29) 

by using the values used in the U-I-GWA output, 
C/u+GWA = W{0) = 3.21 eV and Ju+gwa = 0.9 eV. 
This result of the energy difference is in good agreement 
with the calculated one of about 1.5~2 eV. We should 
note that the gap between majority-spin e^ and minority- 
spin band is not the Hund splitting, but is given by 
the difference in the Coulomb interaction (the Coulomb 



gap); 

-{ii;c(efef )-£;c(ef)} 

= K., + U[J - «^ - Jt, J = W{0) + 1.91JU+GWA 

~ 4.93 eV, (30) 

which can nicely explain the spectrum in Fig. [9l The 
position of the minority spin aig band locates below that 
of the minority spin e^ band which is estimated as 

{£;c(e^Te-Ta-^i)_iJc(efc^T)} 

-{Ec{c-ghl')-Ec{ol^)} 

= + <J - ("U - ^'t^J = WiO) + 0.37 Ju+gwa 
~ 3.54 eV. (31) 

The observed V 3d spectrum in paramagnetic metal- 
lic phase shows a two-peak structurei^ The prominent 
peak at Ep of the width of 0.5 eV in metallic phase dis- 
appears in insulator phase, corresponding to opening the 
gap, and the broader peak with a width of 2 eV does not 
change much among antiferromagnetic and paramagnetic 
insulator phases. The observed V d-bands in insulator 
phases are the lower Hubbard band and may correspond 
to the occupied U+GWA bands. But the present cal- 
culated one has a significantly narrower width of about 
0.7 eV. In fact, the tig single-electron level locates at 
- J^2, = W{0) - 1.17Ju+GWA ^ 2.16 eV below the 
lowest level and would appear between the calculated 
CgCg and O 2p-bands. Therefore, we believe that the 
observed broad occupied V d-bands originate from the 
mixture of e^Bg and Cg states. The screened Coulomb 
interaction can be calculated in GWA and U-I-GWA. 

Exchange interactions in V2O3 can be calculated also 
in the framework of LSDA+U and Jq, = 31 meV (fer- 
romagnetic), J/3 = —29 meV (antiferromagnetic) and 
J^ = 43 meV (ferromagnetic) for the first, second and 
third nearest neighbor pairs, which are consistent with 
observed spin alignment. The spin magnetic moment is 
evaluated to be 1.9 /ie- The vanadium trivalent ion is a 
system with electrons less than half and, therefore, the 
spin-orbit interaction is negative. The resultant orbital 
contribution to the magnetic moment should be negative. 
This is also consistent with experimental observations of 
the total magnetic moment. 

B. Discussion on paramagnetic V2O3 

In the paramagnetic case, the electronic structure cal- 
culated by LDA-I-U method is quite different from that in 
antiferromagnetic case, e.g., metallic and large bonding- 
antibonding splitting of aig orbital which causes the ma- 
jority configuration to be aigCg . With a larger value of U, 
the aig band can move to higher energy region but, even 
if J7 = 4.0 eV, the aig band overlaps still with occupied 
Cg band. 
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The total electron-electron interaction energy can be 
evaluated as |{2K^^ -jt,J + 2<^ +Mt, J = 1.67M^(0)- 
0.67Ju+GWA in the case of paramagnetic metallic config- 
uration (ef )i/3(ef )i/3(e^i)i/3(c-i)i/3(aT^)i/3(aiji/3 

where we assume that two-electrons occupy the non- 
degenerate aig and doubly degenerate e^ bands with spin 
degeneracy without the Hubbard gap. Similarly we can 
get the energy l{{ui^^-j^^^) + ui^^+ut,^} = 1.50^^(0)- 
0.21Ju+GWA in the case of paramagnetic metallic config- 
uration (c^T)i/^(c^T)i/2(g^iji/2(g^i')i/2^ assuming that 

the aig band is lifted up and is left empty. On the 
contrary, if we have the Hubbard gap of e^ bands and 
aig band is empty, then we have a paramagnetic insula- 
tor configuration e^^e^-'- and the total electron-electron 
interaction energy can be estimated as — W{0) — 
0.40Ju+GWA- We can compare the energies of two con- 
figurations where two electrons occupy e^ orbitals, the 
energy difference is (1.51^(0) - 0.21Ju+gwa) - {W{0) - 
0.4Ju+gwa) = 0.5M^(0) + 0.19Ju+gwa- Therefore, the 
latter configuration Sg^Sg^ should be realized. However, 
the GWA and U-I-GWA could not create the Hubbard 
gap in paramagnetic phase and does not give a realis- 
tic feature of electronic structure in paramagnetic phase, 
both insulating and metallic phases. This is a possible 
explanation consistent with LDA-|-DMFTi21 In the an- 
tiferromagnetic case, the situation changes very easily 
because the spin polarization can open the Coulomb gap 
as we have seen. Then the crystal field splitting makes 
the aig band appears in the insulating gap. 

V. SUMMARY 

In summary, we proposed GWA method starting from 
the LSDA+U calculation, named U+GWA, with energy- 



only self-consistent calculation. The on-site Coulomb 
interaction parameter is determined so that the off- 
diagonal elements of the self-energy become small and we 
start GWA with more localized wave functions or a wider 
band gap. We then apply U+GWA to antifcrromagnetic 
NiO and V2O3, where the LSDA wave functions may be 
more extended. The antifcrromagnetic MnO may be a 
system to which GWA can be applied. We have given 
a general criterion for choosing the on-site Coulomb in- 
teraction U and the principles whether we should start 
with LSDA or LSDA-hU. The band gap, W{0) and spec- 
tra for MnO and NiO can be evaluated with an excellent 
agreement with the observed results. On the contrary, 
the spectra of V2O3 is much narrower in U-f-GWA but 
the observed V d-bands may be a mixture of c^c^ and 
single-electron e^ level. 



The method of unique choice of U values has been 
analyzed in details The U values cannot be determined 
within the U-I-GWA method but with the help of analy- 
sis of the XPS and BIS spectra, one can choose a reason- 
able value of U and consistent physical properties can be 
determined which are in excellent agreement with exper- 
imental values. 
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